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The theoretical status of NNLO QCD corrections to the inclusive radia- 
tive B — > X s j decay in the standard model is briefly overviewed. Emphasis 
is put on recent results for three-loop fermionic corrections to matrix ele- 
ments of the most relevant four-quark operators. 

1. Introduction 

The inclusive B — ► X s j decay mode, a flavor-changing-neutral-current 
process and therefore loop-supressed in the standard model (SM), is known 
to be a sensitive probe of new physics. Obviously, deriving constraints on the 
parameter space of physics beyond the SM relies strongly on both accurate 
measurements and precise theory predicitions within the SM. 

Combining measurements of BaBar, Belle and CLEO pQ, the current 
world average for the branching ratio with a cut EUo > 1.6 GeV on the 
photon energy in the i3-meson rest frame reads [2j 

B(B -> ^7)^.6 GeV = (3-55 ± 0.24 +° f ± 0.03) x 10~ 4 , (f ) 

where the first uncertainty corresponds to a combined statistical and system- 
atic error, the second one is due to the theory input in the extrapolation of 
the measured branching ratio to the reference value £^,o, whereas the third 
one is connected to the subtraction of b — > dj events. The overall error 
of the world average amounts to about 7% which is comparable with the 
expected size of next-to-next-to-leading order (NNLO) QCD effects to the 
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perturbative transition b — > X^ arton ^f. Thus, a complete SM calculation at 
this accuracy level is desired. 

To a large extent, the NNLO program has been finished and the latest 
theoretical estimate based on the results [3] 

B(B -» X sl )^ >lsGcV = (3.15 ± 0.23) x 1(T 4 (2) 

is in good agreement with the experimental value eq. (flj). Here, the error 
consists of four types of uncertainties added in quadrature: non-perturbative 
(5%), paramteric (3%), higher-order (3%) and m c -interpolation ambigu- 

ity (3%). 



2. The effective theory framework 

The partonic decay width T(b —* sj) recieves large contributions of log- 
arithms log Myy/ml- Resumming them at each order of a s is most suitably 
done in the framework of an effective low-energy theory with five active 
quarks by integrationg out the top and heavy electroweak fields. The rele- 
vant effective Lagrangian is given by 

£ fr = £qcdxQEd(u, d,s,c,b) + ^£-V t *V tb J2Ci(n)Qi(fj,). (3) 

The usual QED-QCD Lagrangian for the light SM fields is stated in the first 
term whereas the second term gives the local operator product expansion 
(OPE) with Wilson coefficients Cj(//) and operators Qi(fJ>) up to dimen- 
sion six built out of the light fields. Vij denotes elements of the Cabibbo- 
Kobayashi-Maskawa matrix and Gf the Fermi coupling constant. 
The operator basis reads 

Qi,2 = (sTic)(cTib), 

^3,4,5,6 = (sTib)J2 q {Q T iQ), 

Qi = ■^-^rfi b {n){s L a llv b R )F lxv , 

Qs = -^m b ^){s L a^T a b R )Gl u . (4) 

where V and V represent various products of Dirac and color matrices. 
rn b (fj>) is the bottom mass in the MS scheme and the sum runs over all 
light quark flavours q. 

Consistent calculations of T(b — * 57) in the effective framework are per- 
formed in three steps. Teh Wilson coefficients Cj(/Uo) /io ~ My/ are first de- 
termined at the electroweak scale by requiring equality of Green's functions 
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in the effective and full theory at leading order in (external momenta) / M\y. 
Subsequently, the operator mixing under renormalization is computed by 
evolving Cj(/i) from /j,q down to the low-energy scale ^ with help of effec- 
tive theory Renormalization Group Equations (RGE). Finally, the matrix 
elements with single insertions of effective operators are computed. Non- 
perturbative effects appear only as small corrections to the last step, which 
is connected to the heaviness of the bottom quark and the inclusiveness of 
the B — > X s j decay mode. 

As far as the next-to-leading order precision is concerned, this program 
has been completed already a few years ago, thanks to the joint effort of 
many groups (see for eg. [U [5] and references therein). The NNLO calcu- 
lation, which involves hundreds of three-loop on-shell vertex-diagrams and 
thousands of four-loop tadpole-diagrams, is a very complicated task and, as 
already mentioned in the introduction, large parts have already been fin- 
ished. 

Matching the four-quark operators Qi,...,Qg and the dipole operators Qj 
and Qs at the two- and three-loop level, respectively, has been performed in 
[6j[7]- The three-loop renormalization in the {Qi, ■ ■ ■ , Q§} and {Q7, Qs} sec- 
tors was found in [HJ [9], and results for the four- loop mixing of Q\, . . . , Qq 
into Qj and Qs were lately provided in [10] completing the anomalous- 
dimension matrix. The two-loop matrix element of the photonic dipole 
operator Qj was found, together with the corresponding bremsstrahlung, 
in [114 I12j and confirmed in [13j . Moreover, contributions of the dominant 
operators in the so-called large-/3o approximation (O (qq3o)) to the photon 
energy spectrum have been computed in [fi]. Three-loop matrix elements 
of the operators Q\ and Q2 at O (a%Po) and two-loop matrix elements of 
Qi and Qs were found in [15] as expansions in the quark mass ratio rr? c jviv\. 
Recently, we confirmed the findings of [15] on the matrix elements of Qi^ 
and were able to extend the calculation beyond the large-/?o approximation 
by evaluating the full fermionic contributions [IB] . This calculation is briefly 
reviewed below. Furthermore, in [17] , the full matrix elements of Q\ and 
Q2 have been computed in the large m c limit, m c 3> m&, and subsequently 
used to perform an interpolation to the physical region assuming that the 
large-/?o part is a good approximation at m c = 0. This is the source of the 
interpolation ambiguity mentioned beneath eq. ([2]). 



3. NNLO fermionic corrections to the matrix elements of Qi,2 

Matrix elements of Q\ and Q2 constitute a crucial input for the accuracy 
of the current NNLO estimate eq. ([2]). The intention of our recent work |16] 
was the determination of full fermionic corrections to these matrix elements 
to cross-check the results of [15J and, at the same time, to test the valid- 
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Fig. 1. Plots of the 0{a 2 s rif) corrections to matrix elements of Q\ as function 
of z = m\jm\ with fermionic loops of mass M and ji = mi,, (a): M = nib, 
(b): M — m c . The M = case is also shown for comparison. 



ity of the massless approximation used in the large-/?o approximation. Our 
calculation is based on two different techniques, that were applied to the 
master integrals obtained from IBP reduction. In the case of massless quark 
loop insertions into the gluon propagator of the relevant NLO diagrams, all 
integrals have been performed using the Mellin-Barnes (MB) method. The 
MB representations were derived using an automated package [IS] and an- 
alytically continued with help of the MB package [19J. After expanding in 
z = m^./mf ) the resulting coefficients represented as series over residues could 
be resummed with XSummer |20j . In addition, an exact solution through di- 
rect numerical integration keeping the full z-dependence was obtained. This 
procedure was not applicable in the case of massive quark loop insertions 
due to poor convergence. Instead, the method of differential equations as 
a second approach was utilized. Using the fact that the master integrals 
Vi(z,e) (after rescaling by a trivial factor) are functions of e and the mass 
ratio z~ l a system of differential equations has been generated where the 
right-hand side was again expressed through master integrals with the help 
of relations obtained from the reduction. The solution of this system for 
arbitrary values of y proceeded in two steps. First, an expansion in e and y 
for e, y —* was performed and the coefficients were calculated recursively 
up to high powers of y [2T]. Using the resulting high precision values at a 
starting point y <C 1, the unexpanded system was integrated numerically 
up to physical values of y with help of the Fortran package ODEPACK [22] . 
The path was shifted into the complex plane to avoid special points. Fig- 
ure Q] shows the resulting data points together with fits for the renormalized 
matrix elements of Q\ with an internal quark of mass M = 0, m c , m^. In the 
case M = it is evident that the massless approximation overestimates 
the massive result by a large factor. For M = m c this difference is not that 
pronounced but still not negligible. 
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4. Conclusions 

Taking new results for the full fermionic corrections at NNLO into ac- 
count, the branching ratio is enhanced by about one to two percent in 
comparison to [17]. Moreover, an evaluation of bosonic corrections at this 
order, thereby completing three-loop matrix elements, is essential to resolve 
the interpolation ambiguity and to further improve the SM prediction for 
the B — ► X s 7 decay. 
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